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TITLE: Localization of High Speed Vehicles Using Continuous Transmit Waves 
INVENTORS: Evan Frank Berkman and Paul Douglas Koenigs 
Related Applications 

[0001] This application claims priority to, and incorporates by reference, the entire disclosure of 
U.S. Provisional Patent Application No. 60/483,713 filed on June 30, 2003. 

Government Interest 

[0002] The disclosed methods and systems were made at least in part based on function by the 
Naval Undersea Warfare Center under Contract No. N66604-00-D-151 A. The U.S. Government 
may have rights in the disclosed methods and systems. 

Field 

[0003] The systems and methods relate to vehicle detection, and more particularly to vehicle 
localization using continuously transmitted frequency modulated signals. 

Background 

[0004] A sonar (or radar) system that is capable of continuously transmitting and receiving 
Continuous Wave (CW) signals can determine a vehicle's speed by measuring the Doppler shift 
but cannot measure range. If the continuous transmission is frequency modulated and the 
vehicle is stationary, the range to the vehicle may be determined by measuring the frequency 
difference between the transmitted signal and the received echo. This concept is used quite 
extensively in altimeters, police radars, Doppler navigators, and personnel motion detectors 
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(Radar Handbook, McGraw-Hill, M. I. Skolnik, editor). Past techniques to obtain the difference 
frequency have used the transmitted signal waveform as the demodulating signal. 

[0005] However, if the vehicle is moving, the instantaneous difference between the transmitted 
and received echo is dependent on range, range rate and, to a lesser degree, time. The 
interdependency results in range Doppler ambiguity found in conventional continuously 
transmitted frequency modulated (CTFM) sonar and radar systems. While Pulse-Echo (PE) 
sonar/radar systems have been developed for determining speed and range of target vehicles, 
such systems typically require significantly higher peak power than CTFM systems with 
equivalent detection and localization performance because the PE system typically operate at low 
duty cycle whereas as a CTFM system works at or near 100% duty cycle. Thus, CTFM can have 
the advantage of requiring substantially lower transmitted power and correspondingly smaller 
transmitters than equivalent performance PE systems. A CTFM system can also provide faster 
update rates. Thus, a CTFM system that can resolve the range Doppler ambiguity can have 
significant advantages over known PE systems. 

Summary 

[0006] Systems and methods described herein can provide for the selection of multiple, suitably 
designed, demodulation reference signals enabling the measurement of a vehicle's speed without 
ambiguity due to range of separation using continuously transmitted frequency modulated 
(CTFM) signals. The unique solution of vehicle speed can then be used to resolve the range 
Doppler ambiguity found in conventional CTFM sonar and radar systems. This disambiguation 
can enable continuous range and speed estimates of high-speed vehicles with an attendant 
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reduction in transmit power relative to pulse-echo systems with similar measures of performance 
due to the 100% duty cycle of the CTFM system. 

[0007] In one embodiment, a method of localizing a target can include transmitting a continuous 
waveform signal, receiving an echo signal from the target, demodulating the echo signal to a 
baseband, estimating a rangewise Mach number, M n based on combining energy in cells along 
hypothesized tracks of a frequency-time diagram for the echo signal, and determining a first 
range of the target based on the estimated rangewise Mach number. 

[0008] The method can include estimating a rangewise Mach number by choosing a frequency 
detection binwidth, 8f mm , for the cells from the relationship: 



B min is the minimum bandwidth of the continuous waveform signal, and 

Q is the number of cells in a cycle time, T. 

[0009] The method can include estimating a rangewise Mach number using matched filters 
having a minimum baseband frequency resolution, df mim of: 




= B mn = 4M r Q 



, where: 



T 




fo is the initial frequency of the continuous waveform signal, and 



Q is the number of cells in a cycle time, T. 
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[0010] Demodulating can include mixing the echo signal with the continuous waveform signal 
using dual extended duration interlaced demodulation signals or dual frequency band 
demodulation signals. In one aspect, demodulating can include mixing the echo signal with the 
continuous waveform signal using dual frequency band demodulation reference signals to 
demodulate the echo signal to a pass band to obtain a demodulated echo signal, and 
demodulating the demodulated echo signal to bring the demodulated echo signal down to the 
baseband. The demodulated echo signal can be band passed filtered prior to demodulating. 

[0011] Transmitting the continuous waveform signal can include determining a frequency 
separation, transmitting a continuous waveform upsweep signal, and transmitting a concurrent 
continuous waveform downsweep signal, wherein the upsweep and downsweep signals have the 
determined frequency separation. Demodulating the echo signal can include demodulating an 
upsweep echo signal to baseband to obtain an upsweep baseband frequency for the upsweep echo 
signal, and demodulating a downsweep echo signal to baseband to obtain a downsweep baseband 
frequency for the downsweep echo signal. Determining the range can include combining the 
upsweep baseband frequency and the downsweep baseband frequency to determine a range rate 
and a second range of the target. Combining can include determining a sum and a difference of 
the upsweep baseband frequency and the downsweep baseband frequency to determine the range 
rate and second range, respectively, and combining the first and second ranges to localize the 
target. 

[0012] Determining the first range can include estimating a frequency difference, Af, between 
said transmit signal and said echo signal, and solving for the first range using the relationship: 
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c Af 



+ 2M r f 0 +4M r (\-M r )f o T 0 , where 



2/ 0 (l-2M r )L(l-M r ) 



fo is an initial frequency of said continuous waveform signal, and 

T 0 is and observation time less than a cycle time T. 

[0013] A method of localizing a target can include determining a frequency separation for non- 
overlapping frequency reception, transmitting continuous waveform upsweep and downsweep 
signals separated by the determined frequency separation, receiving an upsweep and downsweep 
echo signal from the target, demodulating the upsweep and downsweep echo signals to 
baseband, and combining the upsweep baseband frequency and the downsweep baseband 
frequency to determine a range of the target. Combining the upsweep and downsweep baseband 
frequencies can include determining a sum and a difference of the upsweep baseband frequency 
and the downsweep baseband frequency to determine a range rate and range, respectively, of the 



[0014] The frequency separation can be determined based on maintaining the frequency 

separation at least as large as a potential Doppler shift taken from a relationship 

/J - f 2 - 4M max / 0 , wherein the Doppler shifts of the continuous waveform upsweep and 

downsweep signals are expressed as f u = f x + f 0 t and f d = / 2 + f 0 t , respectively, where M max is 



[0015] The sum of the upsweep baseband frequency and the downsweep baseband frequency can 
be based on the relationship Af d (t) + Af u (t) = -2M r (f 2 + f x ). The difference of the upsweep 



target. 



a maximum rangewise Mach number, and f 0 = 
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baseband frequency and the downsweep baseband frequency can be based on the relationship 

AC(0 " A/ B (0 = +2r 0 / 0 (l - 2M r ) + W r (\ - M r )fy - 2M r (f 2 - f t ) , where r a ^ — , r is the 

c 

range, and c is speed in the medium of concern. 

[0016] Demodulating the upsweep echo signal and the downsweep echo signal can include 
mixing the upsweep echo signal with the continuous waveform upsweep signal, and mixing the 
downsweep echo signal with the continuous waveform downsweep signal, wherein the mixing 
can use dual extended duration interlaced demodulation signals, or dual frequency band 
demodulation signals. 

[0017] A system for localizing a target can include a transmitter to transmit a continuous 
waveform signal, a receiver to receive an echo signal from the target, and a signal processor to 
demodulate the echo signal to baseband, estimate a rangewise Mach number, M r , based on 
combining energy in cells along hypothesized tracks of a frequency-time diagram for the echo 
signal, and determine a range of the target based on the rangewise Mach number. 

[0018] In one embodiment, a system for localizing a target can include means for transmitting a 
continuous waveform signal, means for receiving an echo signal from the target, means for 
demodulating the echo signal to baseband, means for estimating a rangewise Mach number, M r , 
based on combining energy in cells along hypothesized tracks of a frequency-time diagram for 
the echo signal, means for determining a range of the target based on the rangewise Mach 
number. 

[0019] In one embodiment, a system for localizing a target can include means for determining a 
frequency separation, means for transmitting a continuous waveform upsweep signal, /„, means 
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for transmitting a concurrent continuous waveform downsweep signal, fd, said upsweep signal 
and said downsweep signal having said frequency separation, means for receiving an upsweep 
echo signal from said target, means for receiving a downsweep echo signal from said target, 
means for demodulating said upsweep echo signal to baseband to obtain an upsweep baseband 
frequency for said upsweep echo signal, means for demodulating said downsweep echo signal to 
baseband to obtain a downsweep baseband frequency for said downsweep echo signal, and 
means for combining said upsweep baseband frequency and said downsweep baseband 
frequency to determine a range of said target. 

[0020] A computer-readable medium can contain instructions for controlling a computer system 
to localize a target by controlling said computer system to transmit a continuous waveform 
signal, receive an echo signal from the target, demodulate the echo signal to baseband, estimate a 
rangewise Mach number, M n based on combining energy in cells along hypothesized tracks of a 
frequency-time diagram for the echo signal, and determine a range of the target based on the 
rangewise Mach number. 

[0021] The instructions can include instructions to choose a frequency detection binwidth, 5f min , 
for the cells from the relationship: 

fc.-^-^. where: 

B min is the minimum bandwidth of the continuous waveform signal, and 
Q is the number of cells in a cycle time, T. 
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[0022] The instructions can include instructions to estimate a rangewise Mach number using 
matched filters having a minimum baseband frequency resolution, Sf mim of: 



fo is the initial frequency of the continuous waveform signal, and 

Q is the number of cells in a cycle time, T. 

[0023] The instructions to demodulate can include instructions to mix the echo signal with the 
continuous waveform signal using dual extended duration interlaced demodulation signals or 
dual frequency band demodulation signals. In one aspect, the instructions to demodulate can 
include instructions to mix the echo signal with the continuous waveform signal using dual 
frequency band demodulation reference signals to demodulate the echo signal to a pass band to 
obtain a demodulated echo signal, and demodulate the demodulated echo signal to bring the 
demodulated echo signal down to the baseband. The instructions can include instructions to 
band pass the demodulated echo signal prior to demodulating. 

[0024] The instructions to transmit the continuous waveform signal can include instructions to 
determine a frequency separation, transmit a continuous waveform upsweep signal, and transmit 
a concurrent continuous waveform downsweep signal, wherein the upsweep and downsweep 
signals have the determined frequency separation. The instructions to demodulate the echo 
signal can include instructions to demodulate an upsweep echo signal to baseband to obtain an 
upsweep baseband frequency for the upsweep echo signal, and demodulate a downsweep echo 
signal to baseband to obtain a downsweep baseband frequency for the downsweep echo signal. 
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[0025] The instructions to determine the range can include instructions to combine the upsweep 
baseband frequency and the downsweep baseband frequency to determine a range rate and a 
second range of the target. The instructions to combine can include instructions to determine a 
sum and a difference of the upsweep baseband frequency and the downsweep baseband 
frequency to determine the range rate and second range, respectively, and combine the first and 
second ranges to localize the target. 

[0026] The instructions to determine the first range can include instructions to estimate a 
frequency difference, Af, between said transmit signal and said echo signal, and solving for the 
first range using the relationship: 



r = — 



2/ c (l-2M r ) 



4/ +2M r f o +4M r (\-M r )f o T 0 



, where 



_(l-M r ) 

fo is an initial frequency of said continuous waveform signal, and 

To is and observation time less than a cycle time T. 

[0027] In one embodiment, a computer-readable medium can include instructions to control a 
processor to localize a target. The instructions can include instructions to determine a frequency 
separation for non-overlapping frequency reception, transmit continuous waveform upsweep and 
downsweep signals separated by the determined frequency separation, receive an upsweep and 
downsweep echo signal from the target, demodulate the upsweep and downsweep echo signals to 
baseband, and combine the upsweep baseband frequency and the downsweep baseband 
frequency to determine a range of the target. The instructions to combine the upsweep and 
downsweep baseband frequencies can include instructions to determine a sum and a difference of 
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the upsweep baseband frequency and the downsweep baseband frequency to determine a range 
rate and range, respectively, of the target. 

[0028] The instructions to determine a frequency separation can include instructions to maintain 
the frequency separation at least as large as a potential Doppler shift taken from a relationship 
f } -f 2 = 4M max / 0 , wherein the Doppler shifts of the continuous waveform upsweep and 

downsweep signals are expressed as f u = f x + f 0 t and f d = / 2 + f 0 t , respectively, where M max is 



a maximum rangewise Mach number, and f 0 = 



2 j 



[0029] The instructions for determining a sum of the upsweep baseband frequency and the 
downsweep baseband frequency can be based on the relationship 

A/j(0 + A/L(0 = ~~2M r (f 2 + f\) • The instructions for determining a difference of the upsweep 
baseband frequency and the downsweep baseband frequency can be based on the relationship 

A/,(0 - A/ tt (r) = +2r 0 / 0 (l - 2M r ) 4- 8M r (1 - M r )/ 0 / - 2M r (/ 2 - >J) , where r Q = — , r is the 

c 

range, and c is speed in the medium of concern. 

[0030] The instructions to demodulate the upsweep echo signal and the downsweep echo signal 
can include instructions to mix the upsweep echo signal with the continuous waveform upsweep 
signal, and mix the downsweep echo signal with the continuous waveform downsweep signal, 
wherein the mixing can use dual extended duration interlaced demodulation signals, or dual 
frequency band demodulation signals. 
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Brief Description of the Drawings 

[0031] The following figures depict certain illustrative embodiments in which like reference 
numerals refer to like elements. These depicted embodiments are to be understood as illustrative 
and not as limiting in any way. 

[0032] Fig. 1 A illustrates a flow chart of an embodiment of a method of localizing a target 
vehicle using continuously transmitted frequency modulated signals; 

[0033] Fig. IB illustrates a flow chart of an embodiment of a method of localizing a target 
vehicle using upsweep and downsweep continuously transmitted frequency modulated signals; 

[0034] Fig. 2 illustrates a system for continuously transmitting and receiving continuous wave 
signals; 

[0035] Figs. 3A-3C illustrate frequency-time tracks for a transmitted upsweep CTFM waveform, 
the transmitted and echo waveforms, and the demodulated waveform, respectively; 

[0036] Figs. 3D and 3E illustrate frequency-time tracks for a transmitted upsweep CTFM 
waveform with dual extended duration interlaced demodulation reference signals and the 
demodulated waveform, respectively; 

[0037] Figs. 3F and 3G illustrate frequency-time tracks for a transmitted upsweep CTFM 
waveform with dual frequency band demodulation reference signals and the demodulated 
waveform, respectively; 
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[0038] Figs. 3H and 31 illustrate frequency-time tracks for a transmitted upsweep CTFM 
waveform with demodulation to intermediate band with dual frequency band demodulation 
reference signals and the demodulated waveform, respectively; 

[0039] Fig. 4 illustrates a range-time diagram for a moving target; 

[0040] Fig. 5 illustrates a frequency-time track for a transmitted upsweep CTFM waveform and 
the echo waveform for a moving target; 

[0041] Fig. 6 is a range-time diagram illustrating the migration of the echo return start and stop 
times relative to transmission cycle; 

[0042] Fig. 7 illustrates a baseband echo frequency-time track; and 

[0043] Fig. 8 illustrates a frequency-time track for transmitted upsweep and downsweep 
waveforms and the echo waveforms. 

Detailed Description of the Preferred Embodiments 

[0044] To provide an overall understanding, certain illustrative embodiments will now be 
described; however, it will be understood by one of ordinary skill in the art that the systems and 
methods described herein can be adapted and modified to provide systems and methods for other 
suitable applications and that other additions and modifications can be made without departing 
from the scope of the systems and methods described herein. 

[0045] As previously described, current sonar/radar systems capable of continuously 
transmitting and receiving Continuous Wave (CW) signals can determine a vehicle's speed by 
measuring the Doppler shift, but cannot measure range; and Continuous Transmission Wave 
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Frequency Modulated (CTFM) signals can determine range for assumed nearly stationary 
targets. However, current CTFM sonars cannot measure either speed or range unambiguously 
when speed is significant. However, for CTFM signals the time dependent change in the 
difference between transmitted and received frequencies depend only on range rate and this 
dependency can be exploited to determine range rate and then that determination can be used to 
determine range unambiguously. The systems and methods described herein measure the time 
dependent component of the difference frequency by using a family of demodulating signals 
parameterized by range rate covering the expected span of vehicle range rates and either 
searching for maximum output among this family of filters or searching among the filter set for 
the filter that eliminates the time dependency of the frequency difference. The demodulator 
signal parameter, and thus the assumed range rate, that maximizes filter output or produces a 
time independent solution can be used to determine an unambiguous range rate and related 
Doppler shift. The Doppler shift can be subtracted from the difference frequency and the 
resultant can be used to obtain an unambiguous range estimate. 

[0046] A CTFM sonar/radar using a set of speed dependent demodulating signals for tracking 
high speed vehicles can require somewhat more bandwidth (for range resolution) and processing 
channels (for range rate resolution) than a comparable PE sonar/radar. Additionally, design 
considerations for a CTFM sonar/radar system include not overloading the receiver by the 
transmitted signal. However, such considerations can be outweighed by the benefits provided by 
a CTFM sonar/radar system, including lower peak power requirements, a smaller transmit 
section, and faster update rates. To provide a better understanding of the CTFM sonar/radar 
systems and methods described herein, some basic principles of such systems will be described 
further below. For ease of explanation, the systems and methods can be described herein in 
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relation to sonar systems. However, it can be understood that the systems and methods described 
herein can be applicable to other systems that can utilize CTFM signals to localize target 
vehicles or other objects of interest. 

[0047] Continuous transmission (CT) sonar typically can include a cyclical, linear, frequency 
modulated (FM) transmit waveform with a 100% duty cycle. Receivers used in CT sonar 
systems can include a dynamic demodulator and a bank of narrowband filters and narrowband 
energy detectors. The dynamic demodulator can be such that the transmit signal dynamically 
demodulates the received signal to baseband. Thus, as is known, the demodulated echo from a 
target at a constant range is a constant frequency, or tone, with the frequency proportional to 
range. The bank of filters and detectors can correspond to range cells. 

[0048] Waveform and processing parameters for CT sonar can include: 

T : cycle time period, cycle duration, repetition rate for CT; 

B : sweep bandwidth, total bandwidth; 

f 0 : initial frequency; 

Sf : receiver frequency resolution; and 

T 0 : receiver averaging/response time where \/8f <T 0 <T . 

[0049] The transmitted signal can be expressed as: 



*(0 = 2> o C-«r); 



[i] 



n 



s o (t) = mS 0 e m ]*u T (t); 



[2] 



[3] 
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fxmuiO = — 0(0 - fo + fo l '■ instantaneous transmit frequency; and 



In 



[4] 




[5] 



where 




: arbitrary phase; 



\S 0 \ 2 : source level at reference distance; and 



ur(t): unit amplitude window function of time duration T. 

Fig. 3 A can illustrate an example of an upsweep CTFM waveform 300. 

[0050] Consider an ideal case of a stationary target and stationary transmitter and receiver and 
negligible range/time spread in the two-way transmission and in the target or target highlight 
echo response (i.e., "point" target). This ideal may be well approximated for source, receiver, 
and target motion sufficiently slow that the relative range rate between source and target and 
receiver and target produces negligible change in range relative to range resolution during any 
"snapshot" processing interval. In this case, the echo 302 is a time delayed scaled replica of the 
transmitted waveform 300, as indicated in Fig. 3B. The difference in instantaneous frequencies 
between echo and transmitted signals in the latter part of the cycle will be proportional to range: 



0</<r 0 : A/ = (1-^)5 and 



[6] 




[7] 



= fecho'fxmit and 



[8] 
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2r 

t q = — : range information. Where r is range and c is speed in the media. [9] 
c 

[0053] To remove known dynamics of the frequency of the transmitted and echo signals before 
estimating the frequency difference, a dynamic demodulation of the echo to baseband can be 
obtained by mixing with the transmit signal, as in Fig. 3C to obtain demodulated waveform 306. 
For this case, the baseband frequency can equal Af, as shown above. The initial time in each 
cycle that frequency is not proportional to range is referred to in the art as "dead time." Older 
CT systems dealt with dead time by using cycle durations that were much longer than the two- 
way propagation time to the maximum range of operation so that dead time is a relatively small 
fraction of the total cycle time with penalties of more required bandwidth and/or reduced range 
resolution. 

[0054] However, recently developed and documented contemporary methods can eliminate dead 
time at the expense of receiver complexity with no other performance or projector requirement 
penalties. Such methods can include: (1) demodulation to baseband with dual extended duration 
interlaced demodulation reference signals 306, 308 (Fig. 3D), as shown at 310 in Fig. 3E; (2) 
demodulation to baseband with dual frequency band demodulation reference signals 312a, 312b 
(Fig. 3F), as shown at 314 in Fig. 3G; and (3) demodulation to intermediate band with dual 
frequency band demodulation reference signals 316a, 316b (Fig. 3H), as shown at 318 in Fig. 31. 
The preferred methods among these options may depend on total system engineering tradeoffs. 

[0055] In methods (1) and (2), a low pass filter rejects frequencies above frequency B and the 
steady tone is proportional to the range in sum of the two demodulation channels (no dead time). 
For method (2), the first demodulation reference is the transmit Linear Frequency Modulated 
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(LFM) signal and the second demodulation reference is shifted up in frequency by B with regard 
to the first demodulation reference. For both methods (1) and (2), confusing signals can still 
remain within baseband sweep band and bandwidth B, respectively. In method (3), frequencies 
above and below the pass band are rejected. After demodulation to the intermediate band, as 
shown in Fig. 31, a second demodulation brings the echo tonal in the pass band down to the base 
band. 

[0056] For the exemplary stationary target, a number of relationships can be developed. The 

r ■ i • cycle time duration T „ m ^ „ . „ 

number of time resolution bins is: — = = BT . From Equations 7 and 9 

time resolution \/B 

2R r 

and with T = — — , then Af = B . Thus, range resolution , 8r, can be related to 

c R„ 



max 



frequency resolution, Sf, by 5/ = B to yield the number of distinctly resolvable range 



cells as: 



8r 8f 



so that the number of range resolution bins equals the number of frequency resolution bins. 
Noting that Sf^ = — , the maximum number of resolvable range bins within observation time 

T 0 < T can be given by: 



Q m ^=^ L = ^ = BT 0 . [11] 

or- of- 

mm J mm 
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[0057] Conversely, the required frequency resolution, Sf can be driven by the required range 
resolution, 5r, such that Sf = — from Equation 10. In order to maximize the Signal to Noise 

Ratio (SNR), one can select the minimum value of Sf consistent with the cycle duration T, which 

for the ideal stationary target case considered is Sf m[n = ^ . Thus, from Equation 1 1, the 

minimum swept frequency, B mirh is reciprocal to the time resolution needed to obtain the desired 
spatial resolution: B m[n =Q> Sf^ . 



[0058] Considering a minimum frequency resolution for observation time To less than cycle time 

1 ^) T c T 

T, with Sf min = — , then B min = Q • 6/ min = f- = f • — = . — . In other words, if the desired 
T Q T 0 T T 0 2-5r T 0 

response time is less than the cycle time, the minimum total swept frequency B must be 
increased proportional to T/T 0 . Similarly, the swept frequency in the observation interval T 0 is a 
factor 7VTof the total swept frequency. Thus, the swept frequency in a desired response time 



T c 

interval —B must be reciprocal to the time resolution needed to obtain the desired spatial 

T * 2Sr v 

resolution. 



[0059] Thus, for a stationary target and given the maximum range, R max the desired range 
resolution, Sr, and the desired response/update time, one can determine the cycle time, T, the 
sweep bandwidth, B, the initial frequency, fo, and the processing frequency resolution, Sf, from 
the above relationships. As an example, for R max = 6 km, Sr = 10 m, To = 1 sec, and with c being 
1500 m/sec: 



C 1500 
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min rp 
*0 





6000 
10 



= 600; 



B = Q-Sf = 600*1 = 600//z; and 
the frequency sweep rate, BIT = 600/8 = 75//z/sec . 

[0060] As noted above, the above formulations do not account for target motions. It can be 
understood that frequency resolution can be driven by target velocity (Doppler) estimation 
considerations. In addition, one can choose smaller range resolution and hence greater total 
bandwidth to reduce reverberation. The range resolution can be sufficiently small such that the 
target strength of the scattering volume or surface is small compared to the target strength of the 
target of interest. Also, the range resolution can be chosen so as not to over resolve the target or 
highlights or scattering centers on the target of interest. Further, in determining the initial 
frequency, effectively the carrier frequency for B « f 0 , additional considerations beyond those 
known for pulsed sonar can be contemplated. These can include selecting/^ large enough such 
that a sound projector with practical fraction bandwidth can provide the desired absolute 
bandwidth, B, and that demodulation down to intermediate frequency is possible. 

[0061] The above formulations can be extended to a moving target. The movement of the target 
can result in time scale distortion of the echo from the target. With reference to Fig. 4, a plot 400 
of range versus time is shown. The echo 402 received at time r'can be seen to have been 
reflected off a target at a time r earlier. In general, the echo 402 at a time t ' can equal the signal 
404 at a time (/' - 2r) : e(t') = s(t' - 2r) , and cr = r(f' - r) . For first order linearization of the 
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target range/time trajectory at any given time, r(t) -r 0 + v r t , ct = r 0 + v r (t' - r) , and thus 



r _ r 0+ v/ = r(Q 1 
c + v r c 1 + M r ' 



and 



e(t) = s 



t- 



2r(t) 



c(\ + M r ) 



= s 



c 



( 1 > 



[12] 



where M r (t) = r{t)lc = v r (t)/c\ the range wise Mach number. 



[0062] Using the known relationships for LFM transmit signal phase: 



1 I-, 
$xmu (0 = 0) J-\- — Mo 1 = Intfo* + — fot ) 5 instantaneous frequency: 



fxmitiO = ~H0 = f 0 + /<A ^ echo instantaneous phase and frequency can be derived 
In 



from Equation 12 

1 ; 



t~- 



2r(/) 



c(l + M r ) 



2r(/) 



c(l + M f ) 



2r(/) 
c(l + M r ) 



ri-jp 

,1 + M r , 



► ; and 



[13] 



respectively. With r(t) = r 0 + v r t , then Equation 13 can be specialized to: 



f,cho(t) = \fo+fo 



l-M. 



l + M 



r J 



2rjc 
(1 + M r ) 



\-M r 
l + M, 



[14] 



where the terms within the brackets, [ ], represent the second compound Doppler shift on the rate 
of change of frequency and the Doppler modified propagation time delay, respectively. The term 
outside the brackets, { }, represents the first overall Doppler shift on the instantaneous 
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frequency. It can be understood that linearization of the time dependency of r(t) about a current 
target location is valid for sufficiently short observation time, with "sufficiently short" being 
dependent on the actual target kinematics and measurement resolution. 

[0063] From Equation 14, the instantaneous difference frequency for a constant range rate target 
can be determined as: 



4A0 =/<*.« -/_(') = 



J o 



l-M. 



( 



-2M r f 0 - 



4M 



\ + M 



fot 



r J 



1 



(i + M r y 



[15] 



1 B 

where f xmit (t) = — $(t) = f Q + fj and f Q = — , as before. For M r « 1, Equation 15 is well 
2tt T 



approximated by: 



A/(0 = 



- — f 0 (1 - 2M r )-2M r f 0 - 4M r (1 - M r )f 0 t 
c 



(1-Af r ). 



[16a] 



The three terms within the brackets [ ] of Equation 16a respectively represent a range 
proportional term, a Doppler shift term dependent only on range rate, and the LFM variation 
term, which is dependent only on range rate. The first and second time independent terms 
represent an ambiguity between range and range rate. 

[0064] The impact of target motion on the echo frequency can be such that the LFM variation 
term in the base band frequency, Af 9 can cause some loss in processing gain for an exact replica 
correlator, as shown in the closing target example 500 of Fig. 5 for M r < 0. The migration of the 
echo return start and stop times relative to transmission cycle can be illustrated in plot 600 of 
Fig. 6. The echo start/stop times can be derived from the following: 
r 0 +v r (NT + d N ) = cd N ; 

r 0 +v r NT = (c-v r )d N ; 
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_ r 0 +v,NT _ c 
N= (c-v f ) = 




\-M r 



; and 



NT + 2d N = 



r 0 +(\ + M r )NT 
\-M r 



where 2d N is the time of flight for the beginning of the N" 1 cycle transmit waveform to travel to 
the target and return. 

[0065] The ambiguity between the range and range rate in Equation 16a is illustrated with the 
following example: 
fo : 20 kHz; 

target max. range rate : 100 km/hr (55kt) = 28 m/sec M r < 0.02; 
max. range : 6 km; 

sweep bandwidth : 600 Hz and cycle period : 8 sec 
-> dfo/dt = B/T= 600 Hz/8 sec = 75 Hz/sec. 

B 

From Equation 16a and noting that f 0 s-,we can obtain 

B(l-2M r ) ~ 564 Hz; 
2Mf 0 < 800 Hz; and 
4M r (l - M r )B * 48 Hz, 

noting the significant ambiguity between the first range dependent term and the second range 
rate dependent term. 

[0066] It can be seen from Equation 16a that if the LFM variation term, 



fswcep = 4M r (l-M r )f 0 t, 



[16b] 
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can be resolved and accurately measured, then M r can be uniquely determined using Equation 
16a. Knowing M r , Af and the transmit signal parameters, range can be determined 
unambiguously using Equation 16c: 



r - — 



2/ 0 (l-2M,) 



^— + 2MJ 0 + 4M r (1 - M r )f o T 0 



(1-M r ) 



[16c] 



However, in order to avoid spreading the LFM variation term power over more than one 
detection cell at baseband, the minimum bandwidth and/or the time duration of the frequency 
detection bin should be appropriately chosen as described in the following. 

[0067] For observation time To, the sweep of the demodulated echo frequency can be 
approximated by 4M r f 0 T 0 = 4M r B(T 0 /T) = 4M r Q(TjT) • Sf for B = Q • Sf . To determine the 
maximum time the demodulated echo is in a single detection bin, f^ eep I Sf can be set equal to 
one to yield T</T = l/4M r Q. 

[0068] Thus, for high speed targets (large M r ), and/or large numbers of range resolution cells, <5r, 
within the maximum range, R max , the sweep can be over many detection bins. For example, if M r 
= 0.02 and Q = 100, then sweep I Sf = 4M r Q = 8 . Spreading over multiple frequency bins in any 
one time resolution bin can degrade range resolution. Overly wide bins can admit more noise 
than necessary and can be demanding on source bandwidth. 

[0069] Thus, for the application of concern, the baseband echo frequency can sweep over 
multiple detection frequency bins within a continuous transmit cycle for a maximum range rate 
target. This can result in not being able to capture signal energy in cycle time T by integrating 
echo baseband power in a single frequency resolution bin. In order to capture the total signal 

20/558837.4 

-23- 



Atty. Docket No. 03-4039 
3/25/04 

energy in a cycle, one can combine energy in frequency-time bins along range rate hypothesis 
tracks as illustrated in plot 700 of Fig. 7. The number of hypothetical tracks can depend on the 
maximum expected range rate and the range rate resolution. 

[0070] Based on the above constraint, one can state that the detection bin bandwidth can be 
greater than or equal to the sweep rate times the minimum observation time, or 



'4M.B\ f l) _ rTTfB^ 



^5f> MM, 



[17] 



Using the relation B = Q Sf , we can obtain: 



B mio 4M r Q 

mm r-g*- 



= 4M, 



( n\ 



= AM. 



8r 



[18] 



[19] 



Fig. 7 provides an example baseband echo frequency-time track using the above criterion to 
select the frequency and consequent minimum time resolutions so that a maximum range rate 
target does not spread over more than two frequency cells at any one time, as shown by track 702 
and cross-hatched frequency cells 704. 

[0071] Fig. 1 A illustrates a flow chart for a method 100 for localizing a moving target using a 
monostatic CTFM system 200, schematically illustrated in Fig. 2. System 200 can include a 
transmitter 202 for generating a continuous wave signal 204, a directional receiver 206 for 
determining arrival angle and receiving an echo 208 reflected from a target 3, and a signal 
processor 210. The distance 212 between the transmitter 202 and receiver 206 is very small 
compared to the distance 24 between transmitter 202 and target 3 or target 3 and receiver 206. 
The equivalent equations for a bistatic CTFM can be derived in manner completely analogous to 
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that for a Monostatic system. Method 100 can start with the transmission (102) of continuous 
wave signal 204 from transmitter 202 as described by Equation 4. 

[0072] Using the previously described relationships, signal processor 210 can determine (104) 
frequency resolution bins and time resolution bins such that a maximum range rate target does 
not spread over more than two frequency cells at one time and a range resolution can be 
specified using Equation 19. Upon receiving echo 208 as described by Equation 13, as at 106 of 
Fig. 1A, signal processor 210 can demodulate the echo 208 to baseband by mixing echo 208 with 
transmit signal 204 as described by Equation 15, as at 108 of Fig. 1A. The demodulation can 
eliminate dead time using dual extended duration interlaced demodulation signals, or dual 
frequency band demodulation signals. In one embodiment, processor 210 can demodulate echo 
208 to an intermediate frequency, or pass band with dual frequency band demodulation reference 
signals, as indicated in phantom at 108 A. Processor 210 can band pass the desired signal, as at 
108B, and perform a second demodulation 108 that can bring the echo tonal in the pass band 
down to the baseband. 

[0073] Signal processor 210 can estimate, at 1 10, the frequency difference, Af, by measurement 
using standard frequency domain digital signal processing techniques to select the frequency bin 
containing the highest power spectral density. Signal processor 210 can combine energy in the 
cells 704 along hypothesized tracks in the frequency-time diagram (e.g., maximum range rate 
target track 702 and low range rate target track 706, Fig. 7) for echo 208 to obtain an estimate for 
the range wise Mach number, M n using Equation 16b as at 1 12. Combining energy in cells can 
improve signal to noise ratio and improve estimates of M r The target range, shown as 214 or 
216 for a Monostatic case in Fig. 2, can be determined unambiguously using Equation 16c at 1 14 
and method 100 can end. 
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[0074] In one embodiment, a set of LFM matched filters at baseband can capture the total signal 
energy in a cycle. However, it can be noted that the LFM rate of the demodulated echo due to 
target motion can be much less than the LFM rate of transmission. For each initial frequency, or 
range, detection cell, there can be on the order of 2*4M max Q range rate hypothesis tracks or 
matched filters. A very high resolution matched filter can be used to estimate M r from the LFM 
component of the baseband echo frequency. However, since the time-frequency track is roughly 
known from the coarser frequency resolution bins used for detection, as above, there can be a 
manageably small number of long duration, high resolution filters matched to the M r hypotheses 
with the coarser resolution track. Sufficient accuracy in estimating M r can ensure that the error 
in the estimate of the 2M r fo term in Equation 16a is small compared to the selected frequency 
resolution, df\8f» 2s M f 0 , where "»" can be interpreted as "greater than by a factor of two 

or more". Estimating M r by coherent matched filtering of the baseband frequency Af with 
various LFM matched filter hypotheses can have nominal resolution of ~1/(4BT) and, hence, at a 
SNR just sufficient for detection: 



Sm - s Tbt ~* 2£m ' /o = TEf^ 5f>> YEf^ Sf T>> w >>h typically - [20] 



[0075] Since smearing of the echo energy across fixed frequency-time resolution cells is not an 
issue with matched filtering to the LFM component of the baseband echo, the SNR can be 
maximized by selecting the absolute minimum frequency resolution consistent with observation 
time T: 8f = 1/7. However, it can be seen that this selection is not consistent with the above 
criterion for range/Doppler ambiguity resolution. Therefore, a larger value of frequency 
resolution, with attendant loss of SNR, can be used in order to disambiguate range and range rate 
using the matched filters. 
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[0076] From Equation 20, 5f » = ^ , or 

2BT IQ-Sf-T 



5 f » "> 3 \l£= = • [21] 



'2gr in 2\2QT \QT 

Thus, ^ can be chosen as the greater of the value of Sf min from Equation 21, or from the criterion 

AM O 

to prevent spreading of the frequency over multiple bins at one time, Sf^ = — at Equation 

19, as described previously. As an example, given T=5 secs;/ 0 = 20 kHz; and Q = 200, then 

Sf» 3.2 Hz -> 0^= 6.4 Hz. 
Thus, the matched filters described above can provide an alternative method for estimating M r . 

[0077] In one alternative embodiment, concurrent upsweep arid downsweep transmissions can be 
used to resolve the range/range rate ambiguity. The transmissions can be identical except for the 
sign on the rate of change of frequency: 

4£(0 = -To/oO ~ 2M r ) - 2M r f Q - AM r {\-M r )fy , and [22] 
A/,(0 = +r 0 /o(l " 2M r ) - 2M r f 0 + 4M,(1 - M,)/ 0 f . [23] 

The sum and difference of frequencies of the demodulated echoes can disambiguate the range 

and range rate: 

Af d (0 + A/ w (0 = -4M r f 0 , and _ [24] 

46(0 ~ 4f.(0 = +2r 0 / 0 (l - 2M r ) + 8Af r 0 - M,)/ 0 f . [25] 
It can be seen that the sum term of Equation 24 can be independent of range and dependent only 
on range rate. In addition, the difference terms of Equation 25 do not contain a large Doppler 
shift term. Plot 800 of Fig. 8 illustrates the upsweep (802, 804) and downsweep (806, 808) 
transmissions and echoes, respectively. 
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[0078] As described above, the concurrent upsweep and downsweep transmissions method 
doubles the source level and bandwidth requirements, or reduces source level bandwidth 
available for upsweep and downsweep transmission. It can also be noted that since the 
demodulated echo frequencies are time varying, the above summing and differencing operations 
are performed on frequencies measured with time resolution shorter than the full cycle time. 
Further, since the upsweep and downsweep transmissions are in different bands, fo to (fo+B) and 
fo to (fo-B), respectively, the upsweep and downsweep echoes can also be in different bands. 
However, Doppler shifting of the boundary between bands from f 0 to f 0 (l-2M r ) can result in the 
lack of a priori knowledge of the boundary between echo bands, which can prevent separation of 
the echoes from the upsweep and downsweep transmissions by filtering. 

[0079] Separation of the cycle start frequencies by an amount larger than the potential Doppler 
shift (Equation 27) can allow separation of the upsweep and downsweep echoes by filtering: 

fu =f\+U and fd = fi + /o'» where [26] 



f-f 2 =4M max f 0 md f 0 = 



r fx+fi 



[27] 



V 



2 ; 

The echoes can be demodulated to baseband using their respective transmit signals as 
appropriate reference signals. The upsweep and downsweep echo frequencies at baseband can 
be combined to determine range and range rate unambiguously. Substituting from Equation 9 
and Equation 27 into Equations 22 and 23 provides: 

4f.(0 = -— / 0 (1 - 2M r ) - 2M r f - 4M r (l - M r )f 0 t , and [28] 
c 

VA0 = +— /oO ~2M r )- 2M r f 2 + 4M,Q - M r )fy , thus [29] 
c 

¥A0 + A/ u (0 = -2M r (/ 2 + yj) , and [30] 
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4r 



VAO - 4f.(0 = +— /.0 - 2M r ) + 8M,(1 - M r )/ 0 / - 2M r (/ 2 - /j) . 



[31] 



[0080] The frequency resolution for range/range rate ambiguity resolution can be determined 
using an approximation of Equation 25 for an initial analysis: 

¥A0 - A/ u (0 ^ A/ A = +2r 0 f 0 + ZM r f 0 t = 2r 0 / 0 - = 2/ 0 (r 0 - &t) , [32] 



/o 



where Afy = A£,(0 + A/ U (Y) = -4M r f 0 . From Equation 32, we can obtain: 



fo 



2B f 0 



■ + — 



T t 



\2B fo) " V 25 fo j 



Af u , which yields: 



£ r = 



T t) ( T 



[2B /J 



e v& ; and for B«f 0 , 



Noting that with £ r = — £ Tq , then 

<5r 2 Sr 2 8r B Sr B Q 

can confirm that, to a first order of approximation, the range error with a dual upsweep and 
downsweep waveform can be about the same as the desired range resolution. 

[0081] The use of concurrent upsweep and downsweep transmissions can present issues to be 
considered. First, target strength can fluctuate with frequency for realistic time (range) spread 
targets. The transmission loss can also fluctuate with frequency due to multi-path propagation. 
However, the large bandwidth of a CT system can average over frequency fluctuations and can 
provide more robust performance against fluctuating or fading targets than a pulsed CW system. 
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A further consideration can result from the concurrent measurements of the upsweep and 
downsweep echoes. Since the upsweep and downsweep frequencies can be well separated, the 
strength of one of the echoes in a short time frame for a fading or frequency fluctuating target 
can be weak. 

[0082] Fig. IB illustrates a flow chart for a method 150 for localizing a moving target using a 
monostatic CTFM system 200, as schematically illustrated in Fig. 2, with transmitter 202 
concurrently generating upsweep and downsweep continuous wave signals 204. Method 150 can 
start by determining a separation between the cycle start frequencies using Equation 27, as at 
152. Signal processor 210 can determine the separation based on maintaining the separation 
larger than the potential Doppler shift (Equation 27). Transmitter 202 can transmit the 
determined upsweep and downsweep signals using Equation 26, as at 154. It can be understood 
that transmitter 202 can include one or more transmitters for providing separate upsweep and 
downsweep transmission signals. 

[0083] Echoes 208 for the upsweep and downsweep signals (Equation 13) can be received at 
receiver 206, as indicated at 156 of Fig. IB. It can be understood that receiver 206 can include 
one or more receiver processing channels and that separate processing channels can be used for 
down and up Doppler targets, corresponding to closing and opening targets, respectively. Signal 
processor 210 can demodulate (158) the upsweep (Equation 28) and downsweep (Equation 29) 
echoes to baseband using separate appropriate references, as described previously. Signal 
processor 210 can combine the upsweep and downsweep echo frequencies at baseband, as at 
160, to estimate by measurement the sum and difference frequencies. The range can then be 
unambiguously determined at a time, t, using the difference frequency, the range rate and 
Equation 31, as indicated at 162 of Fig. IB, and method 150 can end. 
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[0084] The systems and implementations thereof described herein may not be limited to 
particular components or configurations. For example, signal processor 210 can include one or 
more processors configured to provide the operations described above. In addition transmitter 
202 and receiver 206 can be configured separately or combined together and/or with signal 
processor 210. Separate components can be networked together. System 200 can include 
volatile memory and non-volatile memory for storing data corresponding to input and/or output 
values and other intermediate computations for the methods 100 and 150. Data can be stored in 
one or databases. The memory can include instructions for executing the methods 100 and 150, 
which can be transferred, in the course of operation, from the non-volatile memory to the volatile 
memory and to the processor for execution. 

[0085] Additionally, the flow charts in Figs. 1 A and IB can illustrate an exemplary 
implementation of the features of system 200 and other configurations can be contemplated 
and/or other actions, decisions, and the like can be included. For example, an operator of system 
200 can be given the option to choose method 100 or method 150 for determining range, with 
transmitter 202 providing the appropriate signals and receiver 206 and signal processor 210 
configured for both modes of operation. Additionally, results from method 100 and method 150 
can be compared and/or otherwise combined to assist in improving system 200 accuracy. 

[0086] The systems and methods described herein may not be limited to a particular hardware or 
software configuration; they may find applicability in many computing or processing 
environments. They can be implemented in hardware or software, or a combination thereof. 
Preferably, the methods 100 and 150 can be implemented in computer programs running in one 
or more processors. Implementation languages for the computer programs can include high level 
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procedural languages, object oriented programming languages, assembly and/or machine 
language for communication with a computer system. 

[0087] The computer programs can be stored on a storage media and/or one or more devices 
(e.g., CD-ROM, hard disk, or magnetic disk) that are readable by a general or special purpose 
programmable computer for configuring and operating the computer when the storage medium 
or device is read by the computer to perform the procedures described herein. The system and 
methods can also be considered to be implemented as a computer-readable storage medium, 
configured with a computer program, where the storage medium so configured can cause a 
computer to operate in a specific and predefined manner. 

[0088] References to "microprocessor" and/or "processor" can be understood to include one or 
more microprocessors that can communicate in a stand-alone and/or a distributed 
environment(s), and can thus can be configured to communicate via wired or wireless 
communications with other processors, where such processor can be configured to operate on 
one or more processor-controlled devices that can be similar or different devices. Furthermore, 
references to memory, unless otherwise specified, can include one or more processor-readable 
and accessible memory elements and/or components that can be internal to the processor- 
controlled device, external to the processor-controlled device, and can be accessed via a wired or 
wireless network using a variety of communications protocols, and unless otherwise specified, 
can be arranged to include a combination of external and internal memory devices, where such 
memory can be contiguous and/or partitioned based on the application. Accordingly, references 
to a database can be understood to include one or more memory associations, where such 
references can include commercially available database products (e.g., SQL, Informix, Oracle) 
and also proprietary databases, and may also include other structures for associating memory 
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such as links, queues, graphs, trees, with such structures provided for illustration and not 
limitation. 

[0089] While the method and systems have been disclosed in connection with the preferred 
embodiments shown and described in detail, various modifications and improvements thereon 
will become readily apparent to those skilled in the art. Those with ordinary skill in the art will 
recognize that the arrangement and connectivity of the components shown in the figures are 
merely for illustrative purposes, and can be varied accordingly and components may be 
combined or otherwise reconfigured without departing from the scope of the disclosed systems 
and methods. 

[0090] Accordingly, many additional changes in the details and arrangement of parts, herein 
described and illustrated, can be made by those skilled in the art. It will thus be understood that 
the following claims are not to be limited to the embodiments disclosed herein, can include 
practices otherwise than specifically described, and are to be interpreted as broadly as allowed 
under the law. 
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